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MICHAEL ANSHELEVICH 

Abstract. Free Meixner states are a class of functionals on non-commutative polynomials intro- 
duced in [Ans06]. They are characterized by a resolvent-type form for the generating function of 
their orthogonal polynomials, by a recursion relation for those polynomials, or by a second-order 
non-commutative differential equation satisfied by their free cumulant functional. In this paper, we 
construct an operator model for free Meixner states. By combinatorial methods, we also derive an 
operator model for their free cumulant functionals. This, in turn, allows us to construct a number 
of examples. Many of these examples are shown to be trivial, in the sense of being free products of 
functionals which depend on only a single variable, or rotations of such free products. On the other 
hand, the multinomial distribution is a free Meixner state and is not a product. Neither is a large 
class of tracial free Meixner states which are analogous to the simple quadratic exponential families 
in statistics. 



1. Introduction 

The subject of this paper are states and orthogonal polynomials in non-commuting variables. The 
definition is straightforward. The usual orthogonal polynomials are obtained by starting with a 
measure fi on W^, thinking of M[xi, X2, . . . , Xd\ as a vector space with the (pre-)inner product 

(p,g)= / p(x)g(x)rfMx), 

and applying the Gram-Schmidt procedure to the monomials {xu(i)Xu{2) ■ ■ •a^M(n)}- In the non- 
commutative case, one starts directly with a positive linear functional (state) ip on the algebra of non- 
commutative polynomials M(xi, X2, . . . , a;„), and orthogonalizes the monomials in non-commuting 
variables with respect to the inner product 

(P,Q) = ^[P*(x)Q(x)]. 

Among the general "non-commutative measures" and polynomials orthogonal with respect to them, 
there is a specific class of what is appropriate to call free Meixner states. The classical Meixner class 
[Mei34] consists of familiar distributions — normal, Poisson, gamma, negative binomial, Meixner, 
and binomial — which, somewhat less familiarly, share a number of common properties: their or- 
thogonal polynomials have exponential-form generating functions, they satisfy a quadratic regres- 
sion property [LL60], they generate quadratic natural exponential families [Mor82], they are qua- 
dratic harnesses [Wes93], they are induced by representations of 5u(l, 1) [KVDJ98], and they have 
explicit linearization coefficient formulas [KZOl]. The multivariate Meixner distributions have 
also been investigated, frequently in the guise of quadratic exponential families [Cas96, Pom96], 
though a complete classification is still lacking. Even the infinite-dimensional case was considered 
[SniOO, Lyt03]. 
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In [Ans03], I introduced the free Meixner polynomials, which are a family of orthogonal polynomi- 
als in one variable. The term "free" refers to their relation to free probability, see [VDN92, NS06] 
for an introduction. As a matter of fact, these polynomials have been found independently both be- 
fore and after my work, for example in [Sze22, CT84, Fre98, SYOl, KKN06]. They share a number 
of the Meixner properties listed above, as long they are properly translated into the "free" context, 
see my original paper and also [BB06]. Some of the corresponding distributions also appear in 
random matrix theory, as the limiting distributions in the Gaussian, Wishart, and Jacobi ensembles. 

In [Ans06] I started the investigation of multivariate free Meixner distributions, which are states 
on the algebra of non-commutative polynomials. I continue their study in Section 3. The main 
new tool is to represent these states as joint distributions of certain operators on a Fock space, 
following the more general construction in [AnsOV]. I use this machinery, in combination with 
combinatorial methods, to find explicit formulas for the free cumulants of these states. This provides 
an explanation for the one- variable results in Section 3. 1 of [Ans03] and Proposition 2.2 of [BB06], 
and is the first main result of the paper. The operator representation of the state also allows me to 
handle states that are not necessarily faithful, thus answering a question of the referee of [Ans06], 
where only faithful free Meixner states were considered. 

Having an explicit representation for the cumulants, and being able to handle non-faithful states, 
allows me to describe a number of examples, which is done in Section 4. Among the usual multi- 
variate Meixner distributions, two are familiar, namely the multivariate normal and the multinomial 
distributions. It is well known that the free analog of the multivariate normal distribution is the 
distribution of a free semicircular system, see Section 4.2. The second question treated in this paper 
is: what is the "free" multinomial distribution? I show that the basic multinomial distribution itself 
also belongs to the free Meixner class. In particular, this allows me to calculate the distribution of 
a free sum of d-tuples of orthogonal projections. 

Among states on non-commutative algebras, traces form an important class. The final result in this 
paper provides a way to construct a large family of non-trivial, tracial free Meixner states. These 
turn out to be analogs of simple quadratic exponential families. 

2. Preliminaries 

Variables in this paper will typically come in rf-tuples, which will be denoted using the bold font: 

X = {xi,X2, . . . ,Xd), and the same for z, S, etc. 

2.1. Polynomials. Let M(x) = ]R(xi, X2, .... x^) be all the polynomials with real coefficients in d 
non-commuting variables. Multi-indices are elements u E {1, . . . ,d} for > 0; for |m| = denote 
u by 0. Monomials in non-commuting variables (xi , . . . , a;^) are indexed by such multi-indices: 

Note that our use of the term "multi-index" is different from the usual one, which is more suited for 
indexing monomials in commuting variables. 

For two multi-indices u, v, denote by (m, v) their concatenation. For u with \u\ = k, denote 

{u)"P = {u{k),...,ti{2),ti{l)). 
Define an involution on M(x) via the M-linear extension of 
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A monic polynomial family in x is a family {P^{yi)} indexed by all multi-indices 

oo 

U{uG{l,...,rf}^} 

(with — 1 being understood) such that 

Pjj(x) — xq-\- lower-order terms. 
Note that Pi 7^ P{u)op in general. 
Definition 1. A state on R(x) is a functional 

: R{xi,X2, ...,Xd)^'R 
that is linear, compatible with the *-operation, that is for any P, 

cp[P] = ip [P*] , 
unital, that is if [1] = 1, and positive, that is for any P, 

If [P*P] > 0. 

A state is faithful if in the preceding equation, the equality holds only for P = 0. Unless noted 
otherwise, the states in this paper are not assumed to be faithful. 

The numbers (p [x^] are called the moments of cp. 

A state (f induces the pre-inner product 

{P,Q)^ = ^[P*Q] = {Q,P)^ 

and the seminorm 

l|p|l^ = V^pW^]. 

Throughout the paper, we will typically drop (p from the notation, and denote the inner product and 
norm it induces simply by (•,•), 1 1 • 1 1 . 

We may think of is a "joint distribution" of "random variables" (xi, X2, ■ ■ ■ , Xd). In the remainder 
of the paper, as we did in [Ans06], we will assume that under the state (p, the variables have zero 
mean and identity covariance, 

if [xi\ =0, (fi [x,xj] = 

The last assumption is made primarily so that equation (11) has a clean form. In Section 3.1 we 
briefly describe how to modify the results if that assumption is dropped. 

2.2. Monic orthogonal polynomials states. 

Definition 2. A state has a monic orthogonal polynomial system, or MOPS, if for any multi-index 
u, there is a monic polynomial P^: with leading term x^, such that these polynomials are orthogonal 
with respect to </?, that is, 

(Pa, Pa) = 

for u^v. 

Note that the same abbreviation is used in [DES] to denote a class of multivariate orthogonal poly- 
nomials systems, which is different from ours. 

States that have MOPS were characterized in [AnsOV]. We briefly summarize the results of that 
paper which we will use in the next section. 
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2.2.1. Fock space construction I. Let H = C'^, with the canonical orthonormal basis ei, 62, . . . , e^. 
Define the (algebraic) full Fock space of Ti. to be 



k=0 

Equivalently, J^g,\g{H) is the vector space of non-commutative polynomials in ei, 62, . . . , e^. Fol- 
lowing convention, we will denote the generating vector in = C by O instead of 1. 

For i — 1,2, . . . ,d, define and a~ to be the usual (left) free creation and annihilation operators, 

at (e«(i) <H) e„(2) ® . . . e^^k)) = (g) e„(i) e„(2) . . . ® e„(jfc), 

a~ (e„(i) (g) e„(2) <8) . . . e„(fc)) = (e^, 6^(1) ) e„(2) ... (8) eu(fe). 
For each /c > 2 let C^'^) be an operator 

We think of each C*^*^^ as a c?^ x matrix. Assume that for each k, C^''^ is diagonal and C'^^^ > 0. It 
is convenient to also take C^^^ = I; this corresponds to the identity covariance. Similarly, for each 
i = 1,2, ... ,d and each k >1, let T^^^ be an operator 

Assume that T^^'^^ and C^^^ satisfy a commutation relation (see [AnsOV]). We will denote by % and 
C the operators acting as T^''^ and C^^^ on each component. Finally, let d~ — ajC and 

Xi^at + % + dJ. 

With the appropriate choice of the inner product (•, on the completion TciH) of the quotient of 
J^aigCH), all the operators a^, %, d~ factor through to J^cCH), and each Xj, is a symmetric operator 
on it. 

Theorem 1. (Part of Theorem 2 of[Ans07]) Let (fbea state on ]R(x). The following are equivalent: 

(a) The state (p has a monic orthogonal polynomial system. 

(b) There is a family of polynomials {Pa} such that ip [Pa\ — Ofor allu ^ and they satisfy a 
recursion relation 

Xi = Pi + Bi^^^di, 

d 

•^iPu P(i,u) ~l~ ^ ^ Pi,w,uPw ~l~ ^i,uCui 
w=l 

XiPu — P{i,u) + Bi^ijS^aPiS + (^i,u(l)C'u-P(u(2),u(3),...,u(fe))) 

|«J|=|«| 

with C3 > and, denoting Sj — {s{j), . . . , s{k)), 

k k 
Bi,s,u — Pi,u,s ■ 
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(c) For some choice of the matrices C^^^ and T^^"^ as in Section 2.2.1, the state ip has a Fock 
space representation <fc,{Ti} as 

^ [P(xi, X2, . . . , Xd)] = (O, P{Xi, X^,..., Xa)^) . 

We will also need the following relation between the operators in part (c) and coefficients in part 
(b) of the theorem: 

(1) ^(e„(i) (8) ... (8) eu(fc)) = ^ Bi.s.u.e.^(x) ® ew{k) 

\iS\=k 

and 

(2) C(eu(i) (8) ... (8) e„(fc)) = Cueu{\) <8) . . . <8) e„(jk). 

2.3. Fock space construction II. The following construction is a particular case of the construc- 
tion in Section 2.2.1, but this time we provide full details. As before, let ?Y = C*^, with the canonical 

basis ei, 62, . . . , e^i, denote its (algebraic) full Fock space by JFaig(7i), and the generator of the ze- 
roth component by fi. Let C be an n operator onH^H, which we identify with its x matrix 
in the standard basis. Assume that C is diagonal and 

(3) (/ ® /) + C > 0, 

where / will always denote the identity operator on H. On ^aig('^)> define a new inner product 
using the non-negative kernel 

Kc = (J^^'^-') ® + C)) . . . (/ ® + C) ® J^C^-S)) ((J®2 ^ ^ j®(fc-2)^) 

on each Ti®*^, and denote the completion of jFaig(7Y) with respect to this inner product J^cO^)- If 
the inner product is degenerate, first factor out the subspace of vectors of length zero, and then 
complete. 

For I — 1, 2, . . . , d, let and a~ be the usual (left) free creation and annihilation operators as 
defined in Section 2.2.1. Let Ti,. . . ,Ta be operators on H which we identify with their d x d 
matrices. Assume that each Tj is symmetric and 

{Ti^I)C^C{Ti^I). 

With a slight abuse of notation, we will denote 

Ti^Ti® I®^^-^^ on Ti®^ 

and 

a, = a-(C 7® on 7^®^ 

Note that 

(4) ajn = Tft = = and Si = on TL. 

It follows from the general construction in Section 2.2.1 that all the operators 

Xi — al + a~ +Ti + di 

factor through to J^c(T~(-)- 

Definition 3. The Fock state = v^c,{t,} on ]R(x) determined by such C and Tj is the state 

if [P{xi, X2,..., xd)] = {n, P{x,, X2,..., Xa)n) = {fi, p(Xi, X2,..., Xd)o)^ . 
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2.4. Non-crossing partitions. A partition vr of a set V C Z is a collection of disjoint subsets of V 
(classes of tt), tt = i?2, • • • , B^), whose union equals V . Most of the time we will be interested 
in partitions of {1, 2, . . . , n}. Partitions form a partially ordered set (in fact a lattice) under the 
operation of refinement, so that the largest partition is 1 = ({1, 2, . . . , n}) and the smallest partition 
is = ({1} , {2} , . . . , {n}) . We will use i ~ j to denote that i, j lie in the same class of tt. 

Let NC (V) denote the collection of non-crossing partitions of V, which are partitions tt such that 

i Z i\ j ^ j\ i ^ j,i < j < i' ^ i < f < i'. 

Equivalently, a partition is non-crossing if and only if one of its classes is an interval and the restric- 
tion of the partition to the complement of this class is non-crossing. Non-crossing partitions are a 
sub-lattice of the lattice of all partitions. For each n, let NC{n) denote the lattice of non-crossing 
partitions of the set {1, 2, ... , n}. We will also denote by NCq{V) all non-crossing partitions with 
no singletons (one-element classes), and by NC'iV) all the partitions tt such that 

min V ~ max V. 

Equivalently, partitions in NC\V) have a single outer class — the one that contains both min V and 
max V — in the terminology of [BLS96]. (A class 5 G tt is outer if there do not exist i, i' ^ B, 
j e B with i ~ i' and i < j < i'.) See [NS06] or [Sta97] for more details on the relevant 
combinatorics. 



2.5. Free cumulants. The free cumulant functional R corresponding to a state ip is the linear 
functional on M(x) defined recursively by i? [1] = and for \u\ = n. 



(5) 



■K£NC{n), Be-rv 



.ieB 



which expresses R [x,g] in terms of the joint moments and sums of products of lower-order free 
cumulants. From these, we can form the free cumulant generating function of if via 



(6) 



R{zi, Z2,...,Zd)^^^R [xn] zn, 

n=l \u\=n 



where z— {zi, . . . , za) are non-commuting indeterminates. One can also define R using an implicit 
functional relation involving the moment generating function of cp, see Corollary 16.16 of [NS06]. 

2.6. Words and partitions. In this section, we collect a number of facts that will be useful in the 
proof of the next two theorems. Note that in many places, operators are considered as acting on 
J^aigCH), with a degenerate inner product, rather than on Tc{Ti). 

Lemma 1. Let u be a multi-index indexed by a setV C Z, and W = HiGV (0 ^ word with 
W{i) equal to a'^^yT^^i), a'^^.y or du{i). If 



iev 
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then 

W{mmV) = a-^^,^y^, W{maxV) = a+^^^^), 

(7) yieV,\ {j eV\j> i, W{j) = a-^.^ or WU) = a„o) } | < | { j eV\j> i, WU) = a+ .) } | , 

\{j e = a;(^.) or W^(j) = a„(,)}| = |{j G V\W{j) = a+^.)}| , 

(8) \{j e > = a;(.) oriyO-) = a„o)}| = \{j E V\j > i,W(j) = a^^^}\ 

Proof. This follows from the fact that if e then a+(r;) e j^.j^^^ ^ ^®fc^ ^j^^j 

0^(77), 0^(77) e H^e^-i), and equation (4). □ 

In combinatorics, equation (7) is related to the notion of a Motzkin path. More generally, our 
operator representations are closely related to a common way of representing moments as sums 
over lattice paths [FlaSO, Vie85], but in the multivariate case we find the operator formulation more 
useful. 

Notation 4. Let 

y\;^{u) ^{W ^ W{l)W{2) . . . W{n) satisfying conditions (7) and (8) for y = {1, . . . , n}}, 
and for a general subset y C Z, define {u) similarly. 

Lemma 2. For any multi-index u, partition tt G NCQ{n), tt = (Vi, V2, . . , Vk) and partitions 
(Tj e 7VCo(T^), j — 1,2, . . . ,k, define a word W — /3u(7r; ai, . . . , ak) by 



(9) W{i) 



a^i^j^y i & B & aj,i — maxB, 
VO' i eVj,i^mmVj, 
au(i), i E B E aj,i — min B,i ^ min Vj, 
yTu{i), otherwise. 

Then W e Wn{u), and for each V E tt, W restricted to V is in Wv{u : V), where {u : V) is the 
sub-multi-index ofu indexed by the elements ofV. Moreover, for each u, (5q is a bijection. 

Proof. Let W — Pai'^; ui, . . . , cTfc). Condition (7) for the whole set {1, 2, . . . , n} (respectively, 
for Vj) follows from the definition of /5 and the fact that tt, Ui, . . . , cr^ (respectively, aj) are non- 
crossing. Condition (8) follows from the definition that the minima of the outer classes of the 
partition tt (respectively, aj) are all a~. 

Conversely, let W e W„(m). Let Ac {1,2,..., n}, 

A = {jW(j)7^T„a)}. 

It follows from Proposition 2.13 and Exercise 8.23 of [NS06] that, as long as W restricted to A 
satisfies condition (7), there is a unique non-crossing pair partition tt' e NC{A) such that for any 
B e tt', 

i = minS W{i) = a~^.-^ or du(i), 
i' = max 5 ^ W{i') = aj.,). 
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Note that W{1) = a'^-^y so (1, i') e tt' for some i' > 1. Moreover, W{i' + 1) . . . W{n)Q G H^^, 
so by condition (8), W{i' + 1) = %{i'+iy Thus {i' + 1, j') e tt' for some j', etc. ending with 

(s, n) e tt'. It follows that for any j e {1, . . . , n}, there exist i ^ i' such that i < j < i' and 

= a~^^^. For each j, choose the largest i such that i ^ i', i < j < i', and = a~^^^, 

and require that i ^ j ^ i' . Similarly, for each class K e tt and each j e Vg, choose the largest 

i such that i ~ i' and i < j < i', and require that i ^ j ^ i'. Pictorially, we draw the integers 
1, 2, .... n on a line, draw the pair classes or n' as arcs connecting each i with the corresponding i' 
above the line, and then connect each of the other elements to the arc immediately above it. □ 

Lemma 3. For any W e Wn{u), let (tt; (Ji, . . . , a^) = P^^i^) with tt = (Vi, V2, . . . , Vk). Then 
(O, W{1)W{2) . . . W{n)n) = n ( ^' n ^U)^) ■ 

Proof. Since tt is a non-crossing partition, it has a class V that is an interval, 

V^=[^,^'] = {J|^<J<^'}. 
Since tt restricted to {1, . . . , n} \y is still a non-crossing partition, it suffices to show that 

I i-l n \ / *' \ 

{VL,W{l)W{2)...W{n)n) = ln,\[W{3) W{m\ ln,\[W{mY 

\ i=l i=i'+l / \ 3=i I 

Denote 

T] = W{i' +l)...W{n)Q e H^"^. 
We now show that for any i < j < i', 

W{j) . . . W{n)n ^Cj0v 

for 

Cj^w{j)...w{i')n. 

The proof is by induction. 

W{i')ri = al^^,^r] = eu(i') 0v = {W{i')n) rj. 

If W{j) = a^^.y Cj = eu(j) ® Cj+i- If = Tu{j), then Q = T^j)Cj+i. W{j) cannot equal a'^-y 
Finally, it follows from condition (8) applied toV= [i, i'] that for all j, i < j < i', 

W{j + l)...W{n)Q e H^' 

with s > m. Thus W{j) may equal a„(j) only if s > m -|- 2, otherwise 

w{j)w{j + 1) . . . w{n)n e n^"^. 

But if s > m 2, (j = a^^.)CCj+i. 
It follows that also 

W(i) . . . W{n)n = {W(i) . . . W{i')n) 77 - (O, W(i) . . . W{i')n) 77. 
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Thus 

{n, W{1)W{2) . . . W{n)n) = {n, W{1) ...W{i-1) {Q, W{i) . . . W{i')n) r]) 

= {n, w{i) ...w{i- 1)7]) {n, w{t) . . . w{i')Q) 



I i-l n \ / \ 

^'n^(^) n wm\ui,Y[w{3)n\. 

\ j=l j=i'+l / \ 3=i I 



□ 



Notation 5. For y c Z, denote 

Wy{u) = {We yVy{u)\W{minV) = a-(^,„^),iy(maxV^) = a+_^), 

and none of the other W{i) are equal to a~(j) }• 
The partition tt corresponding to any such W has only one class, tx = {V) E NC{V), and 
/9„-'(>VU^)) = m),a)\a e NC',{V)} - NC',{V). 

Denote 

e{a;V,u) = {Q,f3^{{V),a)Q). 
Lemma 4.IfWe Wn{u) and pz\W) = (vr, a^, . . . , Uk), tt = {V^, I/2, . . . , 14), then 

k I \ k 

(fi, w{i)... w{n)n) = n ( n ^(^)^ ) = n • ^i))' 

j=i \ ieVj I 3=1 

where {u : Vj) is the sub-multi-index ofu indexed by the elements ofVj. 

Proof. This follows from Lemma 3 using Notation 5. □ 

Lemma 5. Suppose that W G W^(m) such that W{1) = a^^^ = aj and W{2) = du{2) = cii- Then 
(5z^{W) = (({1, . . . , n}), a). It follows from condition (9) that for 2 G B G a, we have 2 = mini?. 
Let k — maxS. Then W{k) = a'^f^f.) 



(O, W{1)W{2) . . . W{k) . . . W{n)n) = (n, aJ~a,W{3) . . . a^^,^W{k + l)...W{n- l)a+„)0 
= (n, ajW{k + l)...W{n~ l)a^^^p) (fi, a" 1^(3) . . . a+ . 

Moreover, the map 

{W e W^{u)\W{l) = a-^^^ = aJ, W{2) = a„(2) = di} 



^ U >^{i,fc+i,...,n}((« -.{hk + l,..., n})) X W[,_,y{{u -.{2,..., k})) 

k=3 

n-1 

- U NC',{{1, k + l,...,n})x NC',{{2, ...,k}) 



k=3 

is a bijection. 

Proof. By the same method as in Lemma 3, we deduce that 

W{3) . . .a+,)iy(A:+l) . . .W{n-l)a^^^^n = {W{3) . . . a+.^O) ® (W^(A:+1) . . . M/(n-l)a+^)n) . 
The inner product of this vector with C (e^ ® ej) is the desired expression. □ 



10 M. ANSHELEVICH 

Lemma 6. Suppose that Cij — C{ei® cj) — c for all Let a e NCqIu), 

o = ■ ■ ■ > ^ij(i)} ' ■ ■ ■ ' ■ ■ ■ > ^fej(fc)}) > 

where each class is ordered and bi^i — 1. Then 

k 

e((7; {1, . . . , n} , 1?) = c^~'^ Yl \^u{bi,i),Tu{bi,2) ■ ■ ■ ^«(6i.,(i)-i)e«{6i,^(i)) 



i=l 



Proof. This follows from the definition of P, noting that W{bi^i) = W{1) = a„(i), W{bij) = 
au{bi,i) = ^%{bi,i) j ^ 1' ^iKj(i)) = '^«(6i,^(o)' ^^^^ ° 

3. Main theorems 

Theorem 2. For each i, let 

Si = + Tj + Oj = - 

fee an operator on J-'aigiH). Then the free cumulants of the Fock state (fc,{Ti} fi'om Definition 3 are 
given by the formula R [a; J = 0, 

R[xiP{x)xj] = (ei,P(S)e,-) = (e^, P(S)e,-)^ . 

Proof Since = a^^^^ + + 0^(1), and using Notation 5, for \u\ = n, 

aeNC'o(n) 

= ^ e(cT;{l,...,n},'S:). 

aeNC'oin) 

Similarly, since Xu(i) — a^^) + ^^(i) + a^^^^^ + using Lemma 2 and the preceding equation, 

n k 

= E E E nQ(^-^-(^^^^)) 

fe=l neNCo{n) TjeNC'o(Vj} i=i 
n=iVuV2,...,Vk) j=l,-,k 

= E E n( E Q{'^^■,Vi,{n■.Vi)) 

k=l 7rG7VCo(n) i=l \aieNC'o{Vi) 
n=(Vi,V2,...,Vk) 

n k 

= E E n i^i^-W) ' %^i)(2) ■ ■ ■ %Vi){n-l)e(«:yi)(n) > 

k=l neNCoin) i=l 
7r=iVuV2,...,Vk) 

Thus 

^ [^u] = E n {^{'3-V){l),S(^a:V){2) ■ ■ ■ S(^a:V){n-l)e{u:V){n)) ■ 



FREE MEIXNER STATES 1 1 

Since R [xi] = Lp [xi] = 0, the conclusion of the theorem now follows from the defining relation for 
the free cumulants, namely 



neNC{n) S67r 



n 



.ieB 



□ 



Corollary 3. Each is symmetric and bounded, hence self-adjoint. 

Proof. The symmetry is proved exactly as in Proposition 1 of [AnsOV], or can be deduced from it. To 
prove boundedness, choose m such that ||C|| , UTiH < m. Since |A^C(n)| < 4", and \Q{'K]V,u)\ < 
ml^l, it follows that R [xa] < (4m)l"l and 

Vc,{Ti} [Xu{i)Xu{2) ■ ■ ■^«(n)] < (16m)". 
Thus for each i, \\Xi\\ < 16m. □ 

Notation 6. Let z — (zi, . . . , z^) be non-commuting indeterminates, which commute with x. For 
a non-commutative power series G in z and i = 1, . . . ,d, define the left non-commutative partial 
derivative by a linear extension of = 0, 

Denote by DC = {DiG, . . . , DdG) the left non-commutative gradient. 

For a non-commutative power series G, denote by its inverse with respect to multiplication. 
For a (i-tuple of non-commutative power series G = (Gi, . . . , G^), denote by G^~^^ its inverse with 
respect to composition (which is also a d-tuple). 

Theorem 4. Let if he a state on ]R(x) with a manic orthogonal polynomial system (MOPS), zero 

means and identity covariance. The following are equivalent. 

(a) There exists a non-commutative power series 

F{z) = 1 + {terms of degree > 2) 
and a d-tuple of non-commutative power series U, 

Ui{z) = Zi + higher-order terms, 
such that the polynomials defined via their generating function 

Pa{^>a = F{^) (l - X • U(z)) 

|3|>0 

are a MOPS for ip. 

(b) The polynomials with the generating function 

(10) 5]P5(x)^3= (i-x.(Dit:)<-^)(z) + it:((Di?)<-^)(z)))"' 

|3|>0 

are a MOPS for ip, where R is the free cumulant generating function (6) ofip. 

(c) The free cumulant generating function of (p satisfies, for each i,j, a (non-commutative) 
second-order partial differential equation 

d 

(11) DiDjR(z) = 5ij + J2 B^jDkR(z) + GijDiR(z)DjR(z), 

k=l 
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where Cij > — 1, B^^ — Bl^^, and for each j, k, either B^^ — Ofor all i, or Cju — Ckufar all 
u. 

(d) There is a family of polynomials {Pu} such that tp [P^] = 0/or all u (!) and they satisfy a 
recursion relation 

•^i — Pit 

d 

k=l 
d 

XiP{j,u) = P{i,j,u) + ''^B^jP{k,u) + ^iA^ + Ci,u{l))Pu-, 
k=l 

where C^, Bfj satisfy the same conditions as in part(c). 

(e) There exist symmetric matrices Ti and a diagonal non-negative matrix C with {Ti ® I)C — 
C{Ti ® I) such that ip has a representation <^c,{Ti} as a Fock state of Definition 3. 

We call such states free Meixner states. 

Proof. The equivalence (a) (b) follows from Lemma 4 of [Ans06] and Theorem 3.21 of [Ans04], 
neither of which relied on the assumption that (p is faithful. The equivalence (d) (e) follows from 
the equivalence between the more general Fock space construction and the more general recursion 
relation in Theorem 1. 



(e) =^ (c). By Theorem 2, 



i,i=l ^ ln|>l ^ 



Therefore 



DjR{z) = XI ( {ej, ei) zi + ^ (cj, SqCi) zqZi j 

1=1 ^ |u|>l ^ 



and 



1=1 ^ |3|>1 ^ 

= (cj, Ci) + ^ ( (cj, TiCi) zi+^ (cj, {Ti + ai)S^ei) z^zi j 
1=1 ^ |«|>i ^ 

d / \ d 

= {cj, ei) + {cj, TiCi) zi+^ {ej, TiSaCi) z^zi j + XI 2Z ^^J' ^^^^^i) 

1=1 ^ |u|>l ^ (=1 |u|>l 

d y d d \ 

= {ej, ei) + {ej,Tiek) (e^, e/) zi+^^ {ej, TiCk) (e^, S^ei) z^tzi j 

1=1 ^fe=l |u|>l k=l ^ 



d 

+ ^ 5^ (Cj, Oj-S'ffe;) ZffZi 
l=\ \u\>l 
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where in the last step we have used the fact that {ck} form an orthonormal basis. Using Lemma 5, 
for n> A and u a muki-index on {3, . . . , n — 1} 

{cj^aiSaei) = {n,aJdiSi;tei) 

n-l 

= Yl Yl ^1^) ^2^^) 

*^=3 M/iGWj^ „j (j,(«:{fc+l,...,n-l}),i) 

^2e^{2,...,fc}(*'('^^{3'-.'=})) 

n-l 

where v = {u : {3, . . . ,k — 1})), s = and w — {u : {k + 1, . . . , n})). Thus 

d I d d 

DiDjR{z) = {ej, (^i) + Y \ Y ^^i' + X] X] "^"^'^ ^^^^ 

;=i yk=i \u\>i k=i 

d 

^^iiYY 5Z Y {ej.Wiei) {ei,W2es) z^ZsZijsZi 

d d 

= (cj, Ci) + ^ (cj, T^Cfc) DkR{z) + ^ ^ (ci, ^ire^) (e^-, Sy^ei) z^z^z^zi 

k=l 1=1 u={v,s,w) 

\v\,\iJj\>0 

d 

= i^j, ei) + J2 (^i' ^^e*^) DkR{z) + Ci,DiR{z)D,R{z). 

k=l 

The conditions on the coefficients in part (c) are equivalent to the conditions on the matrices in part 
(e). 

(c) =^ (e). Since the states are assumed to have zero means, the corresponding free cumulant 
generating functions have no linear terms. In that case, a free cumulant generating function R, 
and so the corresponding state ^p, are completely determined by equations (11). Moreover, for any 
choice of {Cij, Bfj] subject to the conditions of part (c), if 

d 
k=l 

and 

C{ei ® Cj) = Cij Ci (g) Cj, 
then those equations are satisfied by i?^^, So the states whose free cumulant generating func- 
tions satisfy the equations in part (c) are exactly the states in part (e). 

(b) =^ (c). (/? has a MOPS, so by Theorem \, (p = V^c,{t,} for some jc'^'^^ ^'^'^^ |- Thus, ip is the joint 
distribution of the operators {Xi, . . . , Xa) on the Hilbert space TciH), with 

Xi^at + % + a7C. 
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Note that since y? has means zero and identity covariance, = and C^^) = I. Using notation 
from Section 2.2.1, and denoting 

{DR)a{2) = i5„(l)i?(z) . . . L>„(|3|)i?(z) 

and 

e-u = e„(i) (8) ... (8) e^d^ri), 

we see that 

(l - X • z + i?(z)) (l) + Y,iDRU7.)ea) 

u 

d 

^n-Y^ZiCi^ R{z)n + Y^{DR)s{z)e^ + i?(z) ^(L'i?)s(z)es 

i=l u S 

^^{DR)u{'z)e^^,u) - J2{^iDiR{z)n + z^{D,R){7){D R)a{-L)e^) 

d d 

-EE ^i{DRU^)%{ea) -EE - /)e„^. 

i u i u 

Since for any function G with zero constant term, 

d 

(12) ^^.AGlzj^Glz), 

i=l 

the preceding expression equals 

d d 

= - E + E(-^-^)"(^)^" "EE 

i=l u i=l u 

d d 

-EE Zi{DR)u{i)%{e^) ^iiDRUz)a-{e - I)e^. 

i u i u 

Using the expansions (1) and (2) from Theorem 1, we now continue the equation as 

d d 

= - E + ^(L'i?)^(z)e^ "EE ^i^^^)u{'^)Hhu) 

i=l u i=l u 

d 

- ^ Zi{Bi^k,jDkR{'z)ej + YBi^(^i,^U),{j,w)DkR{^)^^ 

i,j,k=l u,w 
d 

- ((C(,,,) - l)DiR{z)DjRiz)ej + 5](C(,,,-s) - l)D,R{z)DjR{z)iDRUz)e^j,^ 
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which can be re-organized as 



a a a 

j=l i=l k=l 

d d 

+ E E [DjRiz)iDRUz) -Y^i (SijiDRUz) 

j = l U 1=1 

d 

+ Y.Y.^^AKnU3,n;)DkR{^){DR)a{z) + (C(,,,-„^) - l)D,R{z)D,R{z){DRUz 



k=i aj 

Using equation (12) again, this equals 
(13) 



ZiCj 



= Q + ^ [DiDjR(z) - (^5ij + J2 Bi,k,jDkR(z) + (C(,,,) - l)D,R{z)D^R{z)) 

1,3=1 k=l 
d 

i,j=l u 
d 

+ A,(M),0-,-)^fe^W(^^)-(z) + (C(,,,.„^) - l)DiR{z)D^R{z){DR)a{z) 



k=l w 



If the polynomials {Pa} with the generating function (10) from part (b) are orthogonal, then 

J2 Pu{^){DR)u{7) = (l - X • z + i?(z))"\ 



|u|>0 



and 

(14) 
so that 
(15) 



(l - X • z + i?(z)) + YiDR)u{-2')e^ = Q. 



Equating to zero the coefficient of ziCj in equation (13), we get exactly equation (11) from part (c), 
with = -Bi,A;j and C^j = C{ij) — 1. The conditions on the coefficients follow from the general 
conditions in Theorem 1. 

(e) =^ (b). liip = (pc,{Ti), it follows that in equation (13), B^k,u),{j,w) = BfjSg^a and Qj^it = l+Qj. 
Then that expression equals to 



a a 

[DiDjR(z) - (Sij + J2 B-,DkR(z) + CijDiR(z)DjR(zfj 

' ' ' k=\ 

^3 + ^ ^ 



i,j=l 
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since part (e) ^ (c). So equation (15) holds. As a result, for polynomials with the generating 
function (10), 

|u|>0 |u|>0 

Thus equation (14) holds, and so the polynomials are orthogonal. □ 

3.1. Nontrivial covariance and other extensions. In this section we consider a number of con- 
structions and examples that involve free Meixner states with non-trivial covariances. We still 
assume that they have zero means (for simplicity); if desired, the means pi,...,pd can easily be 
incorporated into the operator model by considering the operators (Xi + pi, . . . , Xd + Pd) instead, 
and the corresponding combinatorics will involve all non-crossing partitions NC (n) rather than the 
non-crossing partitions without singletons NCo{n). 

On the other hand. Theorem 1 requires that for any state with MOPS, the covariance matrices have 
to be diagonal. But now we allow 

(16) iP [x^] = U. 

Note that degenerate variances il^[xf] = are still not permitted. 

3.1.1. Dilations. Let 99 be a Meixner state, and fix positive numbers (ti, ^2, • • • , td). Let ijj be the 
state defined by the IR-linear extension of 

tjj [P{xi, Xd)] = if [P{tlXi, tdXd)] . 

Note that equation (16) holds. It is easy to see that if {Pa} is a MOPS for cp, then 

Qu{Xi, . . . , Xd) = t^PtiiXi/ti., . . . , Xd/td) 

is a MOPS for ip. 

We now briefly state how the results of Theorem 4 get modified for ^p. The generating function 
for the MOPS still has the same "resolvent" form, and any state with MOPS and such a generating 
function arises as a dilation of a free Meixner state. 

[P{xi, . . . , Xd)] = R<p [P{tiXi, . . . , tdXd)] 

which shows how to modify the differential equation satisfied by the free cumulant generating 
function. Similarly, the MOPS satisfy the recursion relation 

d 

XiQ{j,u) — Q{i,j,u) + '^^UBijQ{k,u) + 5ijtl{l + Ci^u(i))Qu-, 

k=l 

where {B^J,Cij} were the corresponding coefficients for cp. Finally, suppose that 99 = ^c,{Ti}, 

represented as the joint distribution of {Xi, X2, . . . , Xd). In the Hilbert space H — with an 
orthonormal basis {cj}, let fi — Uti. On the Fock space Tc{Ti), let 

and extend these definitions C-linearly to , a/, T/ for any f E H. Let 

(17) Xf^ = a+ + Tf^ + al ^ a^C = UX,. 
Then ijj is the joint distribution of (Xj^, Xjj, . . . , XfJ, 

^[p{x)] = {n,p{Xf)n). 
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3.1.2. Free convolution semigroups. Let (/? be a free Meixner state. For i > 0, define a linear 
functional <^^* via its free cumulant functional using relation (5): 

R^mt [P(x)] = ii?^[P(x)]. 
Note that (^™* [xf\ = t. The notation reflects the fact that 

where ffl is the operation of (additive) free convolution; we will not use this property in the paper. 
Using the methods of Theorem 2, it is easy to see that <^^* is a state (and so positive) if any only if 

t + min Cij > 0, 

in other words if t(/ ® /) + C > 0. In particular, by assumption (3), 93^* is always a state for t > 1; 
this is typical behavior for free convolution, as indicated by Corollary 14.13 in [NS06]. (p^ is a 
state for alH > if and only if C > 0; in this case we say that (p is freely infinitely divisible. 

Again, ip^ has a MOPS and the generating function for the MOPS still has the same "resolvent" 
form. 

d 

DiDjR^mt — 5ijt + B^-DkR^mt + {Cij/t)DiR^mt DjR^mt 
k=l 

and 

d 

k=l 

where {-B^j, Cj-,} were the corresponding coefficients for p. Finally, suppose that cp = (pc,{Ti}- On 
the algebraic Fock space jFaig(7i), define an inner product using the kernel 

= f'Kc/t 

on TY®'^, and denote the completion of J^aig{'H) with respect to this inner product J^q' (H). Let 
(18) = at + Ti + ta- + di = a+ + + ta" (/ + C/t) . 

Then (^^* is the joint distribution of (^x[^\ . . . , X^*^) . 

As constructed above, j-'^j^*^! are represented on different Hilbert spaces for different t. We can 

combine this construction with an idea from Section 7.2 of [SniOO] to represent a whole family of 
functionals {</7^|0 < t < l} on a single space. 

A subset 5* C {1. 2. . . . , n — 1} can be identified with an interval partition 7r{S) e Int{n): if 
5 = {i(l),i(2),...,i(A;)},then 

7r= ({l,...,i(l)},{i(l) + l,...,i(2)},...,{i(A;) + l,...,n}). 

Consider the vector space H = H <Si L°°{[0, l\,dx) as a subspace of the Hilbert space H (8) 
L^([0, 1], dx), with the inner product 

(r/ ® /, C ® c/) = iv, I f{x)g{x) dx. 
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On its algebraic Fock space Tax^{H), define the inner product 

((m ® h)® ■■■®{r]i® fi), iCi^gi)^---^ (Cn ® gn))c 



Hn 



J2 (?7l® ■■■®?7n,C'^'(Cl«'---«'Cn)>n ( / 

l,...,n-l> j=l V"^^ 



Sc{l,...,n-1} 
n{S) = (VuV2,...,Vk} 



ieVi 



dx I , 



where S'^ is the complement {1, . . . , n — 1} \S, and 



Complete with respect to this inner product, to get the Hilbert space Tc{H). On this space, define 



operators 












5 /i) ® . 




5/n)) 


= (Cj ® l[o,t)) ® ® /i) (8) . . . ® (r/n ® /n) 




3 /i) «) ■ 




5/n)) 


= (Ci, 77l) ( J fl{x) dx^ (772 (8) /2) ® ■ ■ ■ ® iVn ® /n), 




5 /i) ® . 




§ /n)) 


= (T^r/i ® /ll[0,t)) ® (7^2 ® /2) ® • • • ® (?7n ® /n), 




5 /i) ^ . 


• ® (r/n 







= ((arC(r7i ® 772)) ® (/ll[0,t)/2)) ® (r?3 ® /s) ^ • • • ® (^n ® fn), 

where l[o,t) is the indicator function of the interval [0, t), and let 



By combining Corollary 3 with (a slight modification of) Theorem 6 from [SniOO], it follows that 
each X^^^ is self -adjoint on J^c(H). Note that if all fi — Qi — l[o,t)> then 

iiVl ® /l) . . . (77n ® fn), (Cl ® ^1) ® • • • ® (Cn ® 9n)) c 



5c{l,...,n-l} 

n{3) = {Vi,V2,-yk) 



VnXl 



C/t 



Moreover, each x/*^ restricted to 



is given by the equation (18), and so (^^* is the joint distribution of (^X^^ , . . . , xj^^^ . 
3.1.3. Rotations. Let O — (Oij) be an orthogonal d x d matrix. Let 

/ d d 

O^x = ^ OiiXi, ■ ■ ■ , 5^ OidXi 



and 
(19) 



[P(x)] = [P(O^x)] 



FREE MEIXNER STATES 



19 



We call i{P a rotation of t{P is the joint distribution of (-^/i, . . . , from (17), where we take 

d 

fj = 0{ej) = ^Oijei. 

i=l 

^p^ need not have a MOPS, since the matrix C need not be diagonal in the basis {/i, . . . , fd}- In 
fact, it follows from Lemma 9 of [Ans06] that (^*^ has a MOPS for all O if and only if Cjj = c for 
all i, j, and that in this case ^p^ is also a free Meixner state. It is easy to see that more generally, if 
S d {1, . . . , 0?} and = c for all i, j e S, then c/?*^ is a free Meixner state whenever 0{ek) — 
for all k ^ S. 

3.1.4. Linear transformations. Finally, one can consider a general invertible change of variables 

(d d 
^ ] -^ilXi, . . . , ^ ^ AidXi 
i=l 1=1 

and the corresponding state </7^ defined as in equation (19). (p^ is the joint distribution of operators 

from (17), where we take /, = A{ej) = Yl'i=i ^ij^i- alternative to our definition, one can 

call free Meixner states all states obtained by a linear transformation of a free Meixner state with 
MOPS (compare with [Pom96]). 



4. Examples 

4.1. Free products. In preparation for the examples in this section, for the reader's convenience 
we explain a key notion from free probability. Again, see [VDN92, NS06] for more details. 

Let (fi, . . . ,ipdhe one-dimensional states on M[a;i], . . . , M[a;d], respectively. There is a canonical 
way to define their /ree product state (p on M(xi, . . . , x^). Combinatorially, a natural way to define 

(p is via its MOPS. Let be the MOPS for For a multi-index u, decompose 

i(2) i(k) 

where the consecutive indices vij) ^ v{j + 1), although non-consecutive indices may coincide. 
Then the MOPS {P^r} for ip are defined by 

For example. 

Note that if one considers polynomials in commuting variables and assumes that all v{j) above 
are different, one gets the usual (Cartesian) product of measures. Also, is a free product state 
if any only if the elements xi, 0:2, . . . , are freely independent with respect to 99, in the sense of 
Voiculescu. This can be taken as the definition of free independence; note that for random variables 
independent in the usual probabilistic sense, their joint distribution is a product measure. Finally, 
the crucial property of free cumulant generating functions is their relation to free products: a state 
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(/? is a free product state of </7i, . . . , if any only if the free cumulant generating function of 
decomposes as 

d 

1=1 

This is often stated as the "mixed free cumulants are zero" condition. It is related to the familiar 
property that the Fourier transform of the joint distribution of independent random variables is the 
product of their individual Fourier transforms. 

It is easy to see that free product free Meixner states are exactly the free products of one-dimensional 

free Meixner states, see Remark 6 of [Ans06]. Recall that these one-dimensional distributions, as 
described in that remark. Theorem 4 of [Ans03] and Section 2.2 of [BB06], are known. With 
variance t, they are: the semicircular (free Gaussian) distributions ■^V4:t — x"^ dx, the Marchenko- 

Pastur (free Poisson) distributions + possibly one atom, and more generally 

1 V4(t + c) - {x-hf ^ 

■dx + zero, one, or two atoms, 



27r 1 + {b/t)x + (c/i2)a;2 
depending on the particular values of b, c, t. 

4.2. Semicircular systems. Let C = and all Tj = 0. Then 

Si = at 

and 

for |u| > 1. Thus in distribution, all ~ 0. Only second-order free cumulants of [Xi, X2, . . . , X^) 
are non-zero, and (/? is the distribution of a freely independent semicircular system, the free analog 
of the standard d-dimensional Gaussian distribution. 

4.3. Free Poisson states. Let C = and Ti arbitrary. Then 

Si = at + T, 

and 

R \X{X'[lX (^z; ^n^j) (^^; ^u^j) • 

Thus in distribution, (5*1, 6*2, ... , Sd) ~ (Ti, T2, . . . , T^). It is appropriate to say that in this case, the 
joint distribution cp of (Xi, X2, . . . , Xd) is d-dimensional free Poisson. In [Ans06] we showed that 
if (p is tracial, then (/? is a rotation of a free product of one-dimensional free Poisson distributions. 
Whether or not (p is tracial, the vector Q is cyclic and separating for the von Neumann algebra 

W*{X^,X2,...,Xd). 

4.4. Free product states. For two vectors /, g, denote by Ef^g the corresponding rank one operator, 

Eah)^f{g,h). 

For an orthonormal basis {/«}, Ef.j^ are the corresponding matrix units. For the standard basis 
{cj}, we will denote these simply by E^. In particular, E^ is the orthogonal projection onto e^. 

Let C(ej (g) Cj) = Ci5ij{ei ® ej), and let Tj = biEu. Then Si acts entirely on the subspace 
.7^ci(Span (cj)), on which it equals 

Si = at + bi + Cia~. 
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af + Citti has the centered semicircular distribution with variance Cj (note that on Tc{'H), this 
operator is not symmetric, so its 5 tor-distribution is different from the semicircular one). Therefore 
Si has the semicircular distribution with mean hi and variance q. Also, it follows that 

R ^XiXji^x j\ (cjj Sjj^Cj^ 

unless 

i = u{l) = . . . — u{n) = j. 

In other words, all the mixed free cumulants of {Xi, X2, . . . , X^) are zero. This says precisely that 
their joint distribution (/? is a free product of the distributions of each of Xi, X2, . . . , Xcj. Each of 
these, in turn, is a one-dimensional free Meixner distribution, whose free cumulants are, up to a 
shift of index, the moments of the semicircular distribution with mean 6j and variance Cj. 

4.5. Exponentiated semicircular systems. Let C{ei (g) Cj) = Ci{ei (g) ej) and Ti — bil. Then 

Si = a+ + CittJ + bil, 

again the distribution of 5'; is the semicircular distribution with mean hi and variance q, but now the 
operators Si themselves are freely independent with respect to the state (p, so that their joint distri- 
bution is a free product. The joint distribution of (Xi, X2, . . . , X^) is not a free product (typically, 
not even tracial); it was described in the last section of [Ans06]. 

Note that the preceding two examples make sense for —1 < q < 0, except that one loses the 
interpretation of Si as having a semicircular distribution with variance Cj, and the resulting states 
are not freely infinitely divisible. 

Remark 1. The constructions in the previous two examples coincide in the one-dimensional case. 
That case, and in particular the corresponding free cumulants, were also considered in [Ans03] 
and described completely in [BB06]. Moreover, many one-dimensional free Meixner distributions 
arise as limits in the central and Poisson limit theorems for the ^-transformed free convolution, in 
the sense of [BWOl]; that paper also contains a Fock space construction which coincides with the 
one-dimensional version of the one in Section 2.3. 

4.6. Free multinomial states. It is well known that the Bernoulli distribution 

(1 -p)5o+p5i 

is a Meixner distribution. It was noted in [Ans03] that it is also a (one-dimensional) free Meixner 
distribution. Moreover, the binomial distributions, which are convolution powers 

((1 - p)5o + = E (") (1 - pr v% 

k=0 ^ ^ 

of the Bernoulli distribution, are all Meixner, and the free binomial distributions, which are free 
convolution powers ((1 — p)5o + p5i) of the Bernoulli distribution are free Meixner. In fact, it 

was noted in [BB06] that ((1 — p)So + pSi) are free Meixner for all real t >1. 

It is also well-known that the multinomial distributions are Meixner [Pom96]. In particular, the 
basic multinomial distribution 

(20) Pl^ei + P2Se2 + ■■■ PdSea 

on MJ^ has this property. We now show that it, and so the free semigroup it generates, also induce 
free Meixner states. In this example, it is natural to consider the state with non-trivial means and a 
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non-diagonal covariance matrix; an actual free Meixner state can be obtained from it by an affine 
transformation as in Section 3.1. 

In the Fock space construction of Section 2.3, take dimTi = d — 1 rather than d, and 

C{ei Cj) = — Cj Cj 

for all i.j, so that Cij = —1. In this case the induced inner product on the Fock space jFaig(7i) is 
degenerate, and the vector space factors through to simply Tc{H) = €. ® H. In this example only, 
choose (linearly dependent) vectors {ei\i = 1,2, ... d} vaV. that are not orthonormal, but instead 
satisfy 

{ei,ei) -pi), 

(Cj, Cj) PiPji 

where 

Pi > 0, i = l,2, pi +P2 + ■ ■ ■ +Pd = 1- 

Since these numbers are the covariances of the centered version of the basic multinomial distribution 

(20) , the corresponding matrix is positive semi-definite and so the {ej} can be chosen in this fashion. 

Let 

Ti{ei) = (1 - 2pi)ei, 

and define 

Xi = aj + Tj + a~ + a^C 
as usual, except that = on H. In other words, for Yi = Xi+ pi, 

YiQ = ei+ piQ, 

(21) YiCi = {1 - p^)[e^ + piOl, 

YiCj = -pj[ei + PiQ]. 

Proposition 5. Yi is an orthogonal projection ofC^Ti. onto Span (cj + piQ). These projections are 
orthogonal among themselves and their sum is the identity operator. Their joint distribution with 
respect to the state (pc,{Ti}from Definition 3 is the basic multinomial distribution (20). In particular, 
</? [1^] = Pi. The free cumulant generating function of if satisfies the dijferential equation 

DiDjR = {5ijPi - piPj) + {5ij - Pj)DiR - piDjR - DiR DjR 

= 5ij (DiR + pi)- (DiR + Pi) (DjR + pj). 

Proof. Yi is self-adjoint, its image is Span (cj + Pi^l), and 

Yi(ei +piQ) = (1 -pi)[ei + Pifl] + Pi[ei + PiQ] = + PiQ, 
so it is an orthogonal projection onto its image. 

(ci + piQ, Cj + pi^l) = -piPj + piPj = 0, 

so these subspaces, and therefore projections onto them, are orthogonal. C (BH has dimension d, 
therefore the sum Yli=i is identity. It also follows that 

1 0, otherwise. 
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Therefore, the joint distribution of (Yi, I2, • • • , Yd) with respect to ip is the basic multinomial distri- 
bution. The last part follows from the operator representation. □ 

Definition 7. Free multinomial states are the states {v^^*|t > l}, where if is the basic multinomial 
distribution (20). Note that <^^" is the joint distribution of the sum of n (i-tuples of orthogonal 
projections. 

Remark 2. Since the Bernoulli distribution is both classical and free Meixner, one may conjecture 
that it is in some sense also g-Meixner (for < g < 1, with q — 1 corresponding to the classical 
case and q = corresponding to the free case). The meaning of this term is not well-defined, but see 
for example Section 4.3 of [Ans04]. Indeed, the recursion relation for its orthogonal polynomials 
is of the g-Meixner form 

xPq = Pi+p, 

xP,=P2 + {l-p)Pi+pil-p)Po, 

xPn = Pn+1 + (1 - p)['n]qPn + N«p(l - p){l - [u - l]q)Pn-l 

independently of q, as long as the degree of the polynomial n < 1, which suffices since 

L'{{l-p)5o+p5i) 

is 2-dimensional. One may also hope that its g-cumulant generating function would then satisfy the 
equation 



where Dg is the g-derivative 



and 



l-q)z 

00 ^ 

n=l '■ 

The corresponding recursion for its g-cumulants is 

nl 

Oii+lOCn-i+l 



an+2 = Oin+1 - ^ 

i=0 



(compare with Remark 5.4 of [BB06]), with the initial condition ai = p. Using Maple, it is easy 
to calculate the first 5 cumulants. Unfortunately, the fifth g-cumulant of the Bernoulli distribution 
calculated in this fashion differs from its fifth g-cumulant in the sense of Section 6 of [AnsOl]. 

4.7. Tracial examples. If </? is a state on a non-commutative algebra A, one says that cp is tracial, 
or a trace, if for any x,y e A, 

If [xy] = If [yx] . 

Tracial states play a crucial role, for example, in the theory of von Neumann algebras. 

Lemma 7. Letip = ipc,{Ti} be a free Meixner state, represented as the joint distribution of operators 
(Xi, . . . , Xf). Suppose that ip is tracial. Then for all i, j, 

(22) TiCj = TjCi 
and 

(23) TjTj TjTj^ Cj^Ej^j C^jEj^. 
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SO 



(cj, ei) {ej, Ck) Cm + (e,, TjTkCi) = {ej, e,) (cfc, ei) Cu + (e^, TiTkCj) . 
Using equation (22) and the orthonormality of {e,}, 

(ci, ei) {ej, Ck) Cji + (cj, T^T/efe) = (e^-, e^) (e^, e;) C;^- + (e,, TiTjCk) , 



so 



(e^-, Cfc) Cjiei + Tj-T^Cfc = (e*;, ei) CijCj + T/T,efc 



and 



TjTi — TiTj — CijEji — CjiEij. 



□ 



Example 3. A general theorem of Voiculescu (Proposition 2.5.3 in [VDN92]) implies that the free 
product states from Example 4.4 are tracial. It is also easy to see that any rotation of a tracial state 
is tracial. 

Lemma 8.1f(p is tracial, then ip^^ is tracial for all tfor which it is defined. 

Proof. It is easy to see directly from the defining equation (5) that a state Lp is tracial if any only if 
its free cumulant generating functional is tracial. The combinatorial reason is that if a partition 
TT is non-crossing when points {1, 2, . . . , n} are placed on a line, it is also non-crossing when they 
are placed on a circle. The lemma follows from this fact and the defining relation for c/?™* 



Proposition 6. All free multinomial states are tracial. 

Proof. Since the operators [Yi] defined in equation (21) commute, the basic multinomial distri- 
bution is tracial; in fact, it factors through to a state on commutative polynomials M\xi, . . . , x^] 
corresponding to the basic multinomial measure (20). It follows from Lemma 8 that all the other 



We conclude the paper with three further results on when Meixner states are traces. The first two 
show that under one set of general assumptions, the only tracial Meixner states are the trivial ones, 
namely the rotations of free product states. It generalizes Proposition 11 of [Ans06]. The last one 
provides a way to construct a large class of tracial examples that do not come from free products. It 
generalizes the multinomial example above. 

Proposition 7. Let (/? = v?c,{tj be a tracial free Meixner state with C diagonal as a d x d matrix, 
Cij = 5ijCi. Then (p is a rotation of a free product state. 



□ 



free multinomial states are tracial as well. 



□ 
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Proof. If Cij = SijCi, then equation (23) states that all Tj, Tj commute. Combining this with equa- 
tion (22), we see that moreover, for some orthonormal basis {/i, /2, . . . , fd}, 

d 

(24) T, = J]a,(/„e,)^/,j, 

i=l 

From 

C{Tj I)(ek ei) = {Tj ® J)C(efc ® e^) 

it follows that 

d d 

ai {fi, Cj) {fi, Cfe) {fi, ei) ci{ei (g) e;) = 5^ a, {fi, Cj) {fi, eu) {fi, e^) ci{em e^). 

i=l i,m=l 

Thus for all /c 7^ Z, 



^ tti (/i, e^) (/i, Cfe) {fi, ei) ci = ^Cfc, ai {fi, Cj) Ef.j^ 



ei )ci^ {ek,Tjei) ci = 0. 



It follows that whenever q ^ 0, TjCi e Span (e;), and one can take // = e/. So if 5* = ^ 0}, 
then 

Span {ei\l E S) = Span (/||/ G S") 

is an invariant subspace for all Tj. We can choose an orthogonal transformation O so that 0(ej) = 
fi, in other words 

Cj for i E S, that is q 7^ 0, 
fi for i ^ 5, that is q = 0. 



0(ei) 



Following the comments at the end of Section 3.1.3, the state (p"^, which is the joint distribution of 
{Xf-^^, . . . , XfJ, is still a tracial free Meixner state. From equation (24), 

Finally, C{r] <Si C) — ^ whenever one of rj,( E Span {ei\l ^ S) — Span {fi\l ^ S), so 



C{f ® /,) = 



C{ei (g) Cj) = CiSij if i,j e 5*, 

if one of i,j 5". 



Thus C, {Ti} have the form in Example 4.4, and so ip*^ is a free product state. □ 

Corollary 8. Let (p — (pc,{Ti} be a tracial free Meixner state and Ti — Ofor all i. Then (p is a free 
product state. 

Proof. If Tj = for all i, then it follows from equation (23) that Cij — SijCi, and the preceding 
proposition applies. In this case, a rotation is unnecessary. □ 

Meixner states correspond to quadratic natural exponential families. The following states are free 
versions of simple quadratic natural exponential families in the terminology of [Cas96], where all 
such (classical) families were classified. 
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Proposition 9. Let C be a constant matrix, Cij = c for all i.j. Then the necessary conditions 
(22) and (23) for (p to be a tracialfree Meixner state, namely that {Tj} are symmetric matrices, 

TiCj = TjCi and 

{TiTj - TjTi) = c{Eij - Eji), 

are also sufficient. 

Proof. If c = 0, the result follows from Proposition 7. So we will assume that c 7^ 0. 
For each n, let 

A{n) = {ttG 7VC[)(n)|l -2|. 

For any partition a E NCQ{n)\A{n), define the partition l{a) E A{n) as follows: iflEBEa and 
2 G C G cr, let l{a) be the partition with the same classes as a except that i? U C is a class of /(cr). 
Conversely, any such a can be obtained by starting with n E A{ri) with the (unique) outer class B, 
choosing i E B, 2 < i < n (if it exists), and taking a to have the same classes as tt except that 

Bn{{l}U{t + l,...,n}) 

and 

5n{2,...,z} 

are classes of a. 

For TT G NCo{ri), define the partition p(7r) G NCo{{2, 3, . . . , n, 1}) (where 1 is identified with 
n + 1) by 

J ~ J for J / 1, 
1 ~ J 1 ~ J. 

Clearly 

p(7VC[,H) = {tt G iVCo({2, . . . , n, l})|n ^ 1} , 

and 

= {tt G iVCo({2, . . . , n, i})|2 ^ n 1} . 

In particular, p{A{n)) C A^Co({2, . . . , n, 1}). For any partition a E NC'q{{2, ...,n, l})\p{A{n)), 
define the partition r{a) E p{A{n)) as follows: if 1 E B E a and n G C G a, let r(a) be the 
partition with the same classes as a except that i? U C is a class of r(cr). Conversely, any such a 
can be obtained by starting with tt G p{A{n)) with the (unique) outer class B, choosing i E B, 
2 < i < n (if it exists), and taking a to have the same classes as tt except 

sn({2,...,z}u{i}) 

and 

Bn{i + l,...,n} 

are classes of a. 
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We will use the usual commutator notation [T, T'] = TT - T'T. Using equation (22), 

Tu(2)Ttj(3)Ttj(4) . . . Tu(n-2)Tu(n-l)Gu(n)) 

= {eu 

= (e«(2), [Tu(i),Tu{3)]Tu(i) ■ ■ .Tu(„_i)e„(„)) + . . . 

+ (e„(2), T„(3)Tt,(4) . . . Tu(n-2) , T„(„_i)] eu(n)) 

+ (eu(2), ru(3)r„(4) . . . Tu{n-i)Tu{n)^u{l)) 

Now using equation (23), [Tj, T,] = c{Eij — Eji), 

n-l)e„(n)> 

= c (e„(2), e„(i)) (e„(3), r„(4) . . . T„(„_i)e„(„)) 

+ C (e„(2), r^,(3)e„(i)) (e,,(4), T^,(5) . . . Tu(n-l)eu(n)) 

+ C (e„(2), T„(3) . . . Tu(n-2]Gu{:L)) (<^«(n-l)) (^u(n)) + ■ ■ ■ 

+ (e«(2), ^m(3) • • • Tu(n)eu(l)) 

— c (e„(2), 6^,(3)) (e„(i), T„(4) . . . T„(„_i)e„(„)) 

— C (e^j(2), T„(3)e„(4)) (e„(i),T„(5) . . ■ Tu(n-l)eu{n)) — ■ ■ ■ 

— C (e„(2), T„(3) . . . Tu(n-2)eu(n-l)) {&u{l) , e-u{n)) 

The left-hand-side of the preceding equation is equal to 6(in; {1, . . . , n} , m), where In e NC'Q{n) 
is the partition with a single class. Similarly, using Lemma 6, the preceding equation itself states 
that 

e(i„; {1, . . . , n} , «) = e(p(i„); {2, . . . , n, 1} , p{u)) 

n-2 

+ J]e(({2,...,i,l},{i + l,...,n}) {2,...,n,l},p(ir)) 

i=2 
n-1 

- J]0(({2,...,^},{l,^ + l,...,n});{l,...,n},u), 

i=3 

where for a multi-index u = ((m(1), . . . , ^(n)), we denote by p{u) = (m(2), . . . , ■u(n), «(!)) the 
multi-index on {2, ... , n, 1}. Using the descriptions of /(a), r(cr) at the beginning of the proof, this 
in turn equals to 

= e(p(i„);{2,...,n,i},p(M)) 

+ e>(T;{2,...,n,i},p(ir)) 

(25) r:r(r)=p(i„) 

- e>(^);{i,---,^},^)- 

(T:i((T)=i„ 
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Using Lemma 6 again and equation (25) applied to the class of tt containing 1, we conclude that for 

TT e A{n), 

c-(W-i)e(7r; {l,...,n},u) = c-(W-i)0(p(7r); {2, . . . , n, 1} , p{u)) 

+ c J2 c-(M-'^e(r;{2,...,n,l},p(ir)) 

r:r(r)=p(7r) 

-c Yl c-(H-i)e(<7;{l,...,n},xr), 

fT:/(iT)=7r 

or, since |p(7r)| = |7r| and |r| = |a| = |7r| + 1, 

e(7r; {1, . . . , n} , = e(p(7r); {2, . . . , n, 1} , p{u)) 

+ Yl Q{r;{2,...,n,l},p{u)) 

T:r{T)=p{iT) 

- J2 e((7);{l,...,n},^), 



a:l{a)=Tr 



Therefore 



'S'„(2) . . . Su{ 

= 5^ (e(7r;{l,...,n},«)+ e{a; {1, . . . ,n} ,u)) 



7reA(ri) <T:l{a)=n 

= 5^ (e(p(7r);{2,...,n,i},p(u))+ ^ 0(r; {2, . . . , n, 1} , 

7r£A(n) r:r(r)=p(fr) 

J2 Q{T;{2,...,n,l},p{u)) 

Tep(7VC;,({2,...,n,i})) 
= (eM(2), Su{3) ■ ■ ■ Su{n)^u{l)) ■ 

Using Theorem 2, we conclude that the free cumulant functional R^p, and so <^ itself, is tracial. □ 
Example 4. For d = 2, one can take 



c+1 
1 



To 



For c? = 3, one can take 

/O c 
Ti = c 
\0 0. 

For d — A, one can take 



c 0^ 
= ( c+1 I , 
,0 1, 
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1 oy 



^0 0^ 
1 
,0 1 0. 
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Remark 5. If Cij = c for all the corresponding Fock space is an interacting Fock space in the 
sense of [AB98]. If Cij = c and in addition all Tj = 0, the von Neumann algebras W*{Xi, . . . , X^) 
were described in [Ric06]. 
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